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We have studied biomembranes with grafted polymer chains using a coarse-grained membrane 
simulation, where a meshless membrane model is combined with polymer chains. We focus on 
Q^' the polymer-induced entropic effects on mechanical properties of membranes. The spontaneous 
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curvature and bending rigidity of the membranes increase with increasing polymer density. Our 
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lyj ' simulation results agree with the previous theoretical predictions. 
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I. INTRODUCTION 

Cellular membranes in living cells are a complex and compound system, including various 
kinds of proteins and different lipid compositions. Some of them are "decorated" by sugar 
chains called glycoproteins or glycolipids [l|, |2| . From theoretical and experimental aspects, 
pure lipid membrane in tensionless fluidic state has extensively received much attention for 
a long time [3J]. Pioneer studies on the polymer-membrane compound system, especially 
on membrane with anchored polymers, have been recently reported by several groups both 
theoretically and experimentally jj-l3]. They found that the mechanical properties of lipid 
membrane, such as the bending rigidity k, the saddle-splay modulus k, and the spontaneous 
curvature Cq, can be dramatically changed by the interaction between membrane and the 
surrounding solvent and macromolecules J8|, |9|. Liposomes decorated by grafted polymers 
(PEG etc.) are often used for carriers in drug delivery systems 10|, lUl- Thus, it is important 
to understand the effects of grafted polymers for both fundamental research and medical 
applications. 

In this paper, we explain our simulation model and then briefly describe the effects of 
grafted polymers on membrane properties. We employed a meshless membrane model with 
grafted polymer chains, since it allows large scale simulations and the mechanical properties 
of pure membranes are easy to control 12|, ll3|]. For an isolated anchored polymer chain, 
the membrane properties have been investigated in details by using scaling analysis and 



partition function method 



6|. Here we focus on membranes with a high polymer density. 



where the interactions between polymer chains are not negligible and have significant effects. 



II. MODEL AND METHOD 



In this work, we employ a coarse-grained lipid membrane grafted with a coarse-grained 
flexible linear polymer chain to form a grafted biomembrane system (see Fig. [1]). One 
membrane particle represents a patch of bilayer membrane, and possesses only transnational 
degree of freedom. The particles form a quasi-2D membrane by a curvature potential based 
on moving least-squares (MLS) method 13|. (ii) Polymer particles are linked by a harmonic 
potential, and freely move like a bead-spring behavior with a soft-core repulsion. An end of 
each polymer chain is anchored on one membrane particle with the harmonic potential. 
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FIG. 1: A coarse-grained lipid membrane anchored with a coarse-grained polymer chain to form 
a grafted biomembrane system. The dark-gray particle represent membrane particles. A single 
polymer chain is composed of ten (light-gray) particles which are freely linked as an excluded- 
volume chain. 

We consider a patch of membrane composed of A^mb membrane particles. Among them, 
-^chain membrane particles are anchored by a polymer chain. Each polymer chain consists 
of A^p polymer segments with a connected membrane particle. The membrane and polymer 
particles interact with each other via a potential 

Utot = f/rep + f/mb + Up. (l) 

All particles have a soft-core excluded volume term with a diameter of a. 

f/rep = £ ^ exp [-20(ri,/a - I) + B] Untin./a) (2) 

which Tij is the distance between membrane (or polymer) particles i and j. The diameter a 
is used as the length unit, B = 0.126, and /cut(s) is a C°° cutoff function 
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with n = 12. The factor A in Eq. (3) is determined so that /cut(shaif) = 0.5, which implies 
A = ln(2)(scut/shaif)"' — 1- In Eq- (2), we use the parameters A = 1, and Scut = 1-2. 

A. Meshless membrane model 

The membrane potential t/mb consists of attractive and curvature potentials, 

iVmb 

Umh = ^ X^ f^att(Pi) + KapiiVi), (4) 



where the summation is taken only over the membrane particles. The attractive multibody 
potential is employed to mimic the "hydrophobic" interaction. 

f/att(p^) = 0.25 In {1 + exp [-4(p, - p*)]} - C, (5) 

which is a function of the local density of membrane particles 

Pi = ^/cut(ry/a), (6) 

with n = 12, Shaif = 1-8 [satisfied with /cut(shaif) = 0.5], and Scut = Shaif + 0.3. The constant 
C = 0.25 In [1 + exp(4p*)] is chosen so that f/att = at pj = 0. Here we set p* = 6 in 
order to simulate 2D fluidic membrane. For pi < p*, f/att acts like a pairwise potential with 
f^att = ^2 ^ .^^ /cut(^ij/o"). For Pi > p*, this potential saturates to the constant — C. Thus, 
it is a pairwise potential with a cutoff at higher densities than p*. 

The curvature potential is given by the shape parameter called " aplanarity" , which is 
defined by 

api = ,^j^, (7) 

with the determinant D„ = A1A2A3, the trace T„ = Ai + A2 + A3, and the sum of principal 
minors M^ = A1A2 + A2A3 + A3A1. The aplanarity a^i scales the degree of deviation from 
the planar shape, and Ai, A2, A3 are three eigenvalues of the weighted gyration tensor 

aa^iri) = ^(ftj - OiG)i(3j - (3G)wcv{rij), (8) 

j 

where a,f3 E {x,y,z}. Without loss of generality, we suppose Ai ^ A2 ^ A3. If Ai is the 

minimum eigenvalue representing a relative deviation from the local plane patch formed 

by the neighbor membrane particles around, its corresponding eigenvector is collinear with 

the normal vector n of this plane patch. When the i-th membrane particle has two or less 

neighbor particles within the cutoff distance rcc, they could be localized on a certain plane, 

therefore api = 0. The mass center of local plane patch tq = Xlj '^i'"^cv(^ij)/ Xli ""^cvl^ij), 

where a Gaussian C°° cutoff function is employed to calculate the weight of the gyration 

tensor 
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which is smoothly cut off at Vij = r^c- Here we use the parameters n = 12, rga = 0.5rcc, 
and Tec = 3cr. The bending rigidity and the hne tension of membrane edges are hnearly 
dependent on ka and e, respectively, for k^ > 10, so that they can be independently varied 
by changing k^ and e, respectively. As described in the previous work [IS], the aplanarity 
Opi ~ 9Ai/ (A2 + A3) for the smallest eigenvalue Ai ^ A2,A3. Locally a membrane patch can 
be expressed for small fluctuations in the Monge representation as 

z = zo + -Ci{x-XQf + -C2{y-yoY, (10) 

where x and y coordinates are set to the principal axes. Ci and C2 are the principal 
curvatures of the membrane patch. In the case of cylindrical membrane, Ci = 1/R and 
C2 = 0. By averaging over the local neighborhood with a weight function w{r) where 
r^ = x"^ + y"^, we have 

Ai = {z') - {zf = {Cl + CMt' - {r')f). (11) 

Since Cf + C| = (Ci + Ca)^ - 2C1C2, we easily know 

kaapi ~ k^Xi ~ k^{Ci + C2)'. (12) 

Thus, a linear relation between the bending energy in Helfrich's macroscopic model and our 
mesoscopic parameter is obtained as k oc ka, which is numerically confirmed in the previous 
work |l3|,ll4|. 



The membrane has zero spontaneous curvature without polymers. The details of this 



membrane model are described in Ref. 13j 



B. Polymer Chain 

In each polymer chain, its polymer particles are connected by a harmonic potential. 

Up = /Char ^ Uh^r{ri,i+l), (13) 

chain 

where the summation is taken only between neighboring particles in each polymer chain 
and between the end polymer particles and anchored membrane particles. The harmonic 
potential is given by 

UUr) = l{r-br, (14) 

which 6 is a Kuhn length of polymer chain. We choose b = 1.2a here such that a polymer 
chain stays in the force-free state for rj^j+i = b. 



C. Simulation method 

The NVT ensemble (constant number of particles A^, volume V, and temperature T) is 
used with periodic boundary conditions in a simulation box of dimensions L^ x Ly x L^. 
The dynamics of both membrane and anchored flexible polymers are calculated by using 
underdamped Langevin dynamics. The motions of membrane and polymer particles are 
governed by 

where m is the mass of a particle (membrane or polymer particle) and C, is the friction 
constant. gi(t) is a Gaussian white noise, which obeys the fluctuation-dissipation theorem: 

{9^At)) = 0, (16) 

{9iA't)9jA't')) = '2kBT(6ij6ai3S{t-t'), 

where a, (3 E {x, y, z} and k^T is the thermal energy. We use A^p = 10, £ = 4, ka = 10, and 
/char = 10 through this work. 

III. RESULTS 

A cylindrical membrane, composed of A^mb = 1200, with A'chain anchored polymers outside 
(see Fig. [2]) is used to estimate the polymer- induced spontaneous curvature and bending 
rigidity. For a cylindrical membrane with a radius R and a length L^, the curvature free 
energy is given by 



J[^iC^ + C2-CoY + KC^Q 
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where Ci and C2 are the principal curvatures at each position on the membrane surface, 
and the membrane area A = 27tRLz. The coefficients k, and R, are the bending rigidity and 
saddle-splay modulus, respectively. Biomembranes with lipids symmetrically distributing in 
both leaflets of the bilayer have zero spontaneous curvature Cq = 0. The anchored polymer 
outside induces asymmetry to the membrane. 




FIG. 2: Front and lateral snapshots of a cylindrical membrane with anchored flexible polymers 
outside used for simulations. It contains 1200 membrane particles and 120 polymer chains. 




0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 
a/R 



FIG. 3: Force fz dependence on the radius R of the cylindrical membranes with different polymer 
densities (j) = 1/12, 1/10, and 1/8, respectively. The solid lines are obtained by linear least-squares 
fits. 



The bending energy of the membrane generates a shrinking force along the cylindrical 
axis, 
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FIG. 4: Polymer-induced spontaneous curvature is estimated by the linear fitting method. It 
increases with increasing anchored polymer density. 
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FIG. 5: Bending rigidity is estimated by the linear fitting method. It increases with increasing 
anchored polymer density. 

Since a shorter cylinder membrane has a larger cylindrical radius, it has a lower bending 
energy. Previously, we clarified that this force can be used to estimate the spontaneous 
curvature Cq and bending rigidity k in simulations [ij]. We calculated Co and k for another 
meshless model, which has an orientational degree of freedom and can have a finite Cq. The 
obtained values of Co agree very well with those calculated from an alternative method using 
a curved membrane strip. 

We define a polymer density = A^chain/^mb to quantitatively describe the relation be- 
tween the change of mechanical properties of membranes and the anchored polymer concen- 
tration. Figures [3] shows the axial force fz calculated from the pressure tensor with different 
polymer densities. The force /^ increases linearly with 1/R. Thus, Cq and k of the grafted 



membranes can be estimated from the least squares fits to Eq. (fT8l) . The obtained values of 
Co and k are shown in Figs. H] and O respectively. Both quantities increase with increasing 
anchored polymer density. This quantitatively agrees with the previous theory for a low 
polymer density. When the membrane is curved in the opposite direction from the polymer 
grafting, the polymer chains have more space to move. These entropic effects generate the 
spontaneous curvature of the membranes. The details of our results and comparison with 



the theory will be described elsewhere 15 1 



IV. SUMMARY 

We have studied the effects of grafted polymer on the membrane properties. Anchored 
polymer chains are added to a meshless membrane model. The free energy of polymer 
chains is dominated by the entropic effects induced by the degree of freedom of polymer 
conformation. We confirm that the induced spontaneous curvature and bending rigidity both 
increase as the polymer chain density increases. Anchored polymer chains also affect other 



properties of biomembranes (it will be described in Ref. 15(|). The properties and stability 



of membranes can be controlled by the decoration of membrane with grafted polymers. 
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